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Abstract 

We predict the presence and positions of purely-long-range bound states of 4 He(2s 3 5') + 
4 He(2p 3 P) near the 2s 3 Si + 2p 3 Po.i atomic limits. The results of the full multichannel and 
approximate models are compared, and we assess the sensitivity of the bound states to atomic 
parameters characterizing the potentials. Photoassociation to these purely-long-range molecular 
bound states may improve the knowledge of the scattering length associated with the collisions 
of two ultracold spin-polarized 4 He(2s 3 5) atoms, which is important for studies of Bose-Einstein 
condensates. 
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I. INTRODUCTION 



Photoassociation spectroscopy to purely-long-range molecular bound states has proven 



a powerful tool for examining the ground state wavefunction of a 



kali atoms, and for pre- 



HQS. 



In recent years 



cisely determining s-wave scattering lengths for these systems 
there has been a growing interest in rare-gas atoms, and in particular helium systems. The 
production of Bose-Einstein condensates of spin-polarized metastable 4 He atoms by two dif- 
ferent groups allows one to study the properties of these condensates and provides an 
improved environment for further investigation of cold collisions of metastable 4 He BEL 
Additional valuable information about the metastable 4 He system may be obtainable from 

n 

photoassociation spectroscopy Herschbach et al. |8|] observed bound states that dissociate to 
the 4 He(2s 3 Si) + 4 He(2p 3 P2) atomic limit by detecting ions produced by Penning ionization. 
Broadening of the photoassociation peaks due to Penning ionization did not hinder the ob- 
servation of the vibrational series. There have been no published theoretical or experimental 
studies of the 4 He(2s 3 Si) + 4 He(2p 3 P ,i) purely-lon sr-range bound states. 

In this paper we provide theoretical predictions for binding energies of levels of purely- 
long-range molecular potentials that dissociate to the 4 He(2s 3 Si) + 4 He(2p 3 P 01 ) atomic 
limits. The purely-long-range bound states have inner turning points at internuclear separa- 
tions larger than 150 a (la = 5.29177 x 10~ 11 m) and primarily depend on the dipole-dipole 
C3 coefficient and the fine structure splittings of the 3 P atom. These quantities are all deter- 
mined from atomic properties. The purely-long-range bound states are independent of the 
short-range form of the interaction potentials. Consequently, there may not be a detectable 
ion signal from Penning ionization processes, which occur at small internuclear separations, 
R < 20 ao- We cannot make theoretical predictions for levels that have a short-range contri- 
bution, because the multipole expansion of the potentials is not appropriate for R < 30 a 
and no short-range 4 He(2s 3 S) + 4 He(2p 3 P) potentials exist. For the remainder of the paper, 
we implicitly mean 4 He when referring to helium. 

II. METHOD 

We calculate the molecular eigenfunctions of a multichannel Hamiltonian containing the 
relative kinetic energy operator, twelve non-relativistic Born-Oppenheimer potentials which 
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include retardation effects, the atomic spin-orbit interaction of the 3 P atom, and a term de- 
scribing nuclear rotation. The model used is similar to that applied to the alkali systems . 
The Movre-Pichler model p adapted to 3 S + 3 P systems yields the adiabatic potentials 
obtained from diagonalizing, at each internuclear separation, the multichannel potentials 
that include the dipole-dipole and spin-orbit interactions. In our calculations we have not 
included Penning ionization or spontaneous emission, which can both be represented as an 
optical or imaginary potential. The radial extent of the optical potential for Penning ion- 
ization is localized well inside the inner turning point of the purely-long-range bound states 
and therefore has a minor role to play for these states. Spontaneous emission is present at 
all internuclear separations. However, near the 3 S+ 3 P dissociation limits the vibrational 
level spacing is large compared to the atomic spontaneous emission width of 1.63 MHz, and 
spontaneous emission can therefore be treated perturbatively. 

The multichannel Hamiltonian has been set up as described in Ref. The collision 

between two helium atoms conserves the total molecular angular momentum J and parity 
p. The total angular momentum is given by J = 1 + j, where 1 is the rotational angular 
momentum and j = L+S. Here L and S are the electronic orbital and spin angular momenta 
and A and £ are their respective projections onto the internuclear axis. The projection of j 
onto this axis is Q = A + E. In addition, j = j a + jb, where j a and jb are the total angular 
momenta of the 3 S and 3 P atoms, respectively. Since 4 He has zero nuclear spin, p = + 
and p = — parity states correspond to gerade (g) and ungerade (u) states, respectively. 
The cr = g/u labels refer to the inversion symmetry of the electron wavefunctions about the 
center of charge. Negative-parity bound states have contributions from even partial waves 
I only. Consequently, photoassociation of two ultracold metastable helium atoms, where 
s-wave collisions dominate, can only excite ungerade 3 S + 3 P bound states. 

The twelve Born-Oppenheimer potentials are uniquely labelled by 2S '+ 1 £+ for A = or 
2S+1 U cr for |A| = 1, where the + superscript denotes the symmetry of the spatial electronic 
wavefunction with respect to reflection through a plane containing the internuclear axis. 
For R > 30 a , the dispersion potentials are given by /-u (R/ X)C-ia I R 3 + Cqa/R 5 , where /3A 
is an R- and A-dependent retardation correction |l2j. The quantity A = 3258. 17a , where 
A = A/(27r) and A is the 2s 3 S — 2p 3 P transition wavelength. The C3S coefficient is ±26*3 
and C*3n is ±6*3, where C3 is 6.41022 Eha^ jlj|, Eh is a Hartree, and C3 is proportional to 
the square of the 2s 3 S — 2p 3 P transition dipole matrix element. Following the Physical 
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Reference Database Q|, we assume a 1% uncertainty or "extent of the possible error" for 
this coefficient. The C^\/R 3 contributions to the potentials are attractive for the 
5 S+, 1 1I 9 , 3 I1 M , and 5 U g states, and repulsive for the 3 £+ , 5 S+, 1 Tl u , 3 I1 9 , and 5 U U states. 
The van der Waals coefficients are C 6 s = 2620.76 E h al and C 6 n = 1846.60 E h a^ [l5j. These 



U.UU -Cjh U, 

ies [16|. 



can be calculated from atomic dipole moments and transition frequencies [16j. In turns out 
that the C3/R 3 term is the dominant contribution to the dispersion potential relevant for 
the purely-long-range states. 

The jb quantum number for the 2p 3 P atom, which describes the spin-orbit interaction, 
has the values 0,1, or 2. The 2p 3 Pn and 2p 3 P x levels lie 31.9081 GHz and 2.2912 GHz 



Q. 



above the 2p 3 P 2 level, respectively [lJlll8(. Note that these energy splittings do not satisfy 
the Lande intervals. The rotational Hamiltonian is given by h 2 l 2 /2fiR 2 , where \i = m/2 
is the reduced mass of the system. The atomic mass is m = 4.0026032 u, where 1 u 
= 1.66053873 x 10~ 27 kg. 

We investigate the sensitivity of the bound states to the short-range Born-Oppenheimer 
potentials. No data on the short-range shape currently exists, so we use two different forms 
for the short-range potentials for R < 30 ao- The first of these short-range forms assumes 
a hard wall for all twelve potentials at R = 30 a . For smaller internuclear separations the 
potentials are no longer described by the dispersion potential. For the second short-range 
form, the attractive Born-Oppenheimer potentials are given a more realistic Lennard- Jones- 
like form, with inner turning points between 4ao and 5ao- 

The molecule is best described by the Hund's case (a) coupling scheme at small R, 
where the splitting between the Born-Oppenheimer potentials is large compared to the 
spin-orbit interaction. The molecular states are well described by the quantum numbers J 
and 2S+1 A^r. As the internuclear separation increases and the energy separation between 
the Born-Oppenheimer curves becomes comparable to the atomic spin-orbit interaction, the 
Hund's ) coupling scheme becomes less appropriate, and the independent Hund's 

case (c) coupling scheme is a more suitable choice for describing the long-range states and 
potentials. The Hund's case (c) molecular states are described by the quantum numbers J, 
Q^. For Q = the ± superscript label denotes the symmetry of the electronic wavefunction 
with respect to reflection through a plane containing the internuclear axis. We note that 
for the Q = states, only 0~ and 0+ exist for even J, while the 0+ and 0~ states exist for 
odd J. The Hund's case (c) states are a linear combination of the Hund's case (a) states. 
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Table |l] gives the Hund's ) states that contribute to a given Hund's case (c) state. 

The purely-long-range states discussed in this paper are in the radial regime described 
by the Hund's case (c) coupling scheme. We will discuss various approximations to the 
full multichannel Hamiltonian used in calculating these states. The adiabatic Movre-Pichler 
model described previously labels each purely-long-range state by Q^. An improvement on 
this model is the multichannel Hamiltonian that includes the Born-Oppenheimer potentials 
and the spin-orbit interaction. In addition, we discuss two ways of including the nuclear 
rotation. These approximate Hamiltonians are useful for assigning vibrational quantum 
numbers to the bound states calculated using the full multichannel model, and the adiabatic 
Movre-Pichler model is useful for confirming the assignment of Hund's case (c) labels to the 
bound states. 

The adiabatic Movre-Pichler potentials are shown in Figs. and El They are obtained 
using an extension of the model defined in Ref. jjj. This extended Movre-Pichler model is 
adapted to the 3 S + 3 P system and, in addition to the usual dipole-dipole and the spin-orbit 
interactions, includes the Cq/R 6 dispersion term and retardation. As mentioned previously, 
the fine structure splittings for the He(2p 3 P) atom do not satisfy the Lande intervals, 
which results in a coupling between molecular singlet and quintet states. For large R, the 
potentials in Figs. ^ and 121 dissociate to the three 3 P fine structure limits. For smaller R the 
potentials correlate to attractive and repulsive Born-Oppenheimer potentials, as indicated 
in the figures. We find a number of purely-long-range potential wells in the adiabatized 
potentials dissociating to the 2s 3 S\ + 2p 3 Po and 2s 3 Si + 2p 3 P\ atomic limits that can 
support bound states. In Figs.[Uand|2]the l g and 0+ states dissociating to the 2s 3 Si + 2p 3 P 
limit provide an example of such wells, while the wells in the states dissociating to the 
2s 3 5i + 2p 3 P 1 limit are too shallow to be observed on the scale of the present figure. There 
are no purely-long-range states dissociating to the 2s 3 Si + 2p 3 P 2 limit. 



Dy multichannel bound state calcu- 
3]. Long-range bound states are not 



The eigenstates of the Hamiltonian are determinec _ 
lations in a manner similar to that described in Ref. 
"true" bound states of the system when their energy lies above the 2s 3 Si+2p 3 P 2 dissociation 
limit. They are bound states embedded in a continuum and consequently can predissociate. 
In our method a hard wall is placed at a large internuclear separation, R > 5000 ao- Con- 
sequently, the calculated states are either bound states or "box states" , which characterize 
the continuum. These two types of states are distinguished by visual inspection of a multi- 
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channel eigenvector \ty(R)) = Y,i^i(R)\i) or the square root of the density (Y,i l^iiR)] 2 ) 1 ^ 2 , 
where the sum is over a finite number of molecular spin states The function ipi is the 
radial component for spin state which is characterized by the total spin, parity, and other 
angular momenta. The two types of states can also be distinguished by their sensitivity to 
the location of the hard wall boundary condition. The predissociation of the bound states 
will be further discussed in the next section. 



III. RESULTS 



Four sets of purely-long-range ro- vibrational states are found in the multichannel analysis 
and the binding energies for these levels are given in Table |H] We report here only levels 
with binding energies greater than 1 MHz even though states bound by less than twice the 
atomic line width, 3.25 MHz, may not be clearly observable. The levels are characterized 
using Hund's case (c) labels. The l g and 0+ potentials dissociating to the 2s 3 Si + 2p 3 P 
limit and the 0~ and 2 U potentials dissociating to the 2s 3 Si + 2p 3 P 1 limit have purely-long- 
range states. Each of these potentials has no more than 6 vibrational levels. These levels are 
labelled by v and exhibit a rotational progression when J > Q. In addition, 0+ (0~) states 
only exist for J odd (even). The uncertainty in the binding energy due to changing the short- 
range Born-Oppenheimer potentials is negligible. The last significant digit of the binding 
energy indicates the predissociation width, an estimate of which is obtained by comparing 
the energies for varying locations of the hard wall at large R. We find the predissociation 
width is on the order of, or smaller than, the natural line width. Gerade/ungerade symmetry 
states can only be formed by photoassociation via odd/even partial wave collisions. At the 
internuclear separations of relevance to the states described here, retardation corrections 
that introduce dipole-forbidden coupling between gerade and ungerade states [3] are small. 

As shown in Table ITT1 for the lowest allowed rotational level of each state, the 1% extent 
of the error on C3 produces a change of less than 1 MHz for the l g bound states dissociating 
to the 2s 3 5i + 2p 3 P limit and the 0~ and 2 U bound states dissociating to the 2s 3 5i + 2p 3 Pi 
limit. For the 0+ bound states dissociating to the 2s 3 5*i + 2p 3 P limit, changes of up to 4 
MHz are found for the v = to v = 3 levels. The higher vibrational levels are shifted by less 
than 1 MHz. These changes are obtained using a fixed location of the hard wall at large R. 
It is necessary to include retardation as it leads to a small correction in the binding energy. 
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For the deepest purely-long-range potential, the 0+ dissociating to the 2s 3 5'i + 2p 3 P limit, 
the retardation correction lowers the bound state energy of the J = 1, v = ro- vibrational 
level by approximately 7 MHz. The complete Cq/R 6 contribution to the bottom of this 0+ 
well at R « 190 ao is approximately 0.4 MHz. Therefore the shift in the binding energy due 
to any uncertainty in Cq is small compared with the natural line width. 

We look at three examples to clarify the ro-vibrational and Hund's case (c) assignments 
of the levels reported in Table ITT1 In Fig. El we show the binding energies for J = 2, v = to 
v = 3 ro-vibrational levels of the 2 U state dissociating to the He(2 3 Si)+He(2 3 Pi) limit for 
four approximations and the exact multichannel calculation. The four approximations are: 
(A) the adiabatic extended Movre-Pichler approximation, (B) a non-rotating multichannel 
calculation including the Born-Oppenheimer potentials and the atomic spin-orbit interac- 
tion, (C) an adiabatic calculation of the full J = 2 Hamiltonian where the Coriolis coupling 
between different fl^ states has been set to zero, and (D) a multichannel calculation of the 
full J = 2 Hamiltonian where the Coriolis coupling between different states has been set 
to zero. The binding energies labelled by (E) are obtained using a multichannel calculation of 
the full J = 2 Hamiltonian. All four approximations predict the correct binding energies to 
within 10%. Experiments are sufficiently accurate that the full multichannel model must be 
used. For the v = level, the 4.9 MHz difference between approximations (A) and (B) is due 
to the positive diagonal non-adiabatic coupling, which depends on the Born-Oppenheimer 

n 

and spin-orbit interactions 20]. The 21 MHz discrepancy between the two adiabatic approx- 
imations (A) and (C) for the v = level can qualitatively be understood by the perturbative 
inclusion of the h 2 l 2 / (2/ii? 2 ) rotational Hamiltonian. For this Hund's case (c) 2 U state the 
l 2 of the rotational Hamiltonian simplifies to J(J + 1) — 2Q 2 + j(j + 1) = J(J + 1) — 2, 
since j = Q = 2 near the He(2 3 Si)-|-He(2 3 Pi) dissociation limit. Clearly, for the lowest ro- 
tational level J = 2, a non-zero rotational correction occurs. A non-pertubative calculation 
is needed for a quantitative understanding of the discrepancy between the two adiabatic ap- 
proximations (A) and (C) as rotational distortions are important. Similar to the difference 
between approximations (A) and (B), the difference between approximations (C) and (D) is 
also due to a positive diagonal non-adiabatic coupling, which now depends on the rotational 
Hamiltonian as well as the Born-Oppenheimer and spin-orbit interactions. In this case, the 
comparison of approximation (D) with (E) shows that the coupling between Q states shifts 
the J = 2, 2 U levels by no more than 0.5 MHz. Similar comparisons can be obtained for all 
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of the purely-long-range states. 

In Fig. |3] we show the square root of the density for exact J = 2 multichannel wave- 
functions assigned with the 2 U and 0~ symmetries. Also shown are the 2 U and 0~ adiabatic 
potentials dissociating to the 2s 3 S\ + 2p 3 Pi limit, as obtained using the approximation (C) 
described previously. These bound levels are purely-long-range because the square roots of 
the densities show no amplitude for R < 250 do- In addition, the square roots of the densities 
calculated using the full multichannel and adiabatic models do not deviate appreciably. A 
clear vibrational progression from the v = to the v = 1 level of the 2 U state is visible. 

The square root of the density calculated using the full Hamiltonian for three J = 1, l g 
bound levels and two J = 1, l g adiabatic potentials dissociating to the 2s 3 Si + 2p 3 P\ limit 
are shown in Fig. The two potentials are calculated using approximation (C). The "deep" 
l g potential has a double-well structure with a minimum at « 315 Oq. The other "shallow" 
potential has a minimum at ~ 530 a . These levels have not been tabulated in Table HP 
because they are sensitive to the short-range potentials. We present them here as they 
are seemingly purely-long-range except for the penetration of the wavefunction through 
the barrier of the double-well potentials. The peak of the barrier lies slightly above the 
2s 3 Si+2p 3 Pi dissociation limit. The "shallow" potential is not a purely-long-range potential 
because of significant non-adiabatic mixing with the "deep" l g potential. This short-range 
dependence leads to an uncertainty in the level locations and possible broadening by Penning 
ionization processes. In addition, the square root of the density for the v = level has an 
oscillatory behavior for R > 700 a®. The periodicity of these oscillations corresponds to the 
periodicity of the square root of the density for box states in close proximity to the v = 
level. For these l g states the amplitude of the oscillations is primarily determined by the 
rotational coupling. These oscillations indicate significant predissociation. 

IV. CONCLUSION 

We have predicted the binding energies of all four sets of purely-long-range ro- vibrational 
levels of 4 He(2s 3 S) + 4 He(2p 3 P) potentials dissociating to the 2s 3 Si + 2p 3 Pq,i fine structure 
limits. A comparison of the results calculated using the full multichannel and approximate 
models demonstrates that the full multichannel Hamiltonian is required for accurate com- 
parison with experiment. Retardation and the van der Waals Cq interaction are found to give 
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rise to small corrections in the binding energies. We show that a measurement of the binding 
energy to at least 1 MHz leads to a 1% determination of the dipole-dipole C3 coefficient. 

The current measurements of the 5 S+ scattering length of colliding metastable 4 He atoms 
are 300 ± 150 a [5] and 380 ± 190 ao |4j. This implies that a node in the scattering wave- 
function can be found between 150 ao and 600 ao [3]. The inner and outer turning points 
of our purely-long-range bound states approximately fall within this range of internuclear 
separations. Consequently the intensity patterns as a function of vibrational level in pho- 
toassociation spectroscopy may provide an improved measurement of the scattering length. 

The experimental and theoretical binding energies for the 0+ state dissociating to the 



2s 3 Si + 2p 3 Pi limit reported in Ref. 



21] came to our attention upon the completion of this 



paper. Our theoretical predictions are within the experimental uncertainty for the v — 1 — 4 
levels and for the v — level our theoretical prediction of —1418.1 MHz is 2 MHz greater 
than the upper bound on the experimental binding energy —1430 ± 10 MHz. 
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TABLE II: Purely-long-range bound state energies in units of MHz, relative to the dissociation 
limit of the Hund's case (c) state, v and J are the vibrational and rotational quantum numbers. 
The bound states for J < 3 are calculated, and the blank columns indicate that the states do not 
exist for particular J. The last column, AE V , indicates the change in binding energy of the lowest 
rotational level for a +1% change in the C3 coefficient. All binding energies are calculated using 
the full multichannel model. 



Symmetry^dissociation limit 


V 


J = 


J = 1 


J = 2 


J = 3 


AE V 


l g ^2s 3 S 1 + 2p 3 P 







-207.66 


-170.92 


-115.73 


0.62 




1 




-43.80 


-28.52 


-7.97 


0.47 
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-4.12 






0.13 


0+ ^2s 3 S 1 + 2p 3 P 







-1418.1 




-1212.7 


2.4 
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2 




-253.12 




-174.5 


2.5 




3 




-79.65 




-41.6 


1.35 




4 




-18.30 




-4.6 


0.53 




5 




-2.59 






0.13 


0- ^2s 3 S 1 + 2p 3 P 1 





-18.27 




-7.77 




0.08 


2 u ^2s 3 S l + 2p 3 P l 









-191 


-167 


0.5 




1 






-72 


-57 


0.6 




2 






-21.5 


-14.4 


0.2 




3 






-4.7 


-2.2 


0.1 
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FIG. 1: Adiabatic extended Movre-Pichler potentials of gerade symmetry as a function of inter- 
nuclear separation. The zero of energy is the He(2s 3 Si)+He(2p 3 P2) atomic limit. The potentials 
are labelled by Cl^. At short-range we have indicated with ellipses the 25+1 A+ correlations. 




FIG. 2: Adiabatic extended Movre-Pichler potentials of ungerade symmetry as a function of inter- 
nuclear separation. The zero of energy is the He(2s 3 Si)+He(2p 3 P2) atomic limit. The potentials 
are labelled by fi„ . At short-range we have indicated with ellipses the 25+1 A+ correlations. 
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FIG. 3: The binding; energy of the v — to v — 3 vibrational levels of the purely-long-range 2 U 
potential dissociating to the He(2s 3 S'i)-|-He(2p 3 Pi) limit. Various approximations to the complete 
multichannel model are shown. Working from a basic model and increasing in complexity to the 
complete multichannel calculation: (A) the adiabatic extended Movre-Pichler approximation, (B) 
a multichannel calculation including the Born-Oppenheimer potentials and the atomic spin-orbit 
interaction, (C) an adiabatic calculation of the full J = 2 Hamiltonian where the Coriolis coupling 
between different f2„ states has been set to zero, (D) a multichannel calculation of the full J = 2 
Hamiltonian where the Coriolis coupling between different fij states has been set to zero, and (E) 
the complete J = 2 multichannel calculation. 
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FIG. 4: The square root of the density for the J = 2, v = and v = 1 ro-vibrational levels of the 
2 U state and the J = 2, v' = ro-vibrational level of the 0~ state. Note that for the lower two 
levels the square root of the density curves have been inverted to maintain the distinguishability of 
the curves. The solid grey curves represent the adiabatic J = 2, 2 U and 0~ potentials dissociating 
to the He(2s 3 Si)+He(2p 3 Pi) limit. The binding energies of the levels are illustrated by dashed 
horizontal lines. 
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FIG. 5: The square root of the density for the J = 1, v = and v = 1 ro- vibrational levels 
of the deep l g state and the J = 1, v' = ro- vibrational level of the shallow l g state. Note 
that for the lower two levels the square root of the density curves have been inverted to maintain 
the distinguishability of the curves. The solid grey curves represent the adiabatic J = 1, l g 
potentials dissociating to the He(2s 3 Si)+He(2p 3 Pi) limit. The density shows significant amplitude 
for internuclear separations less than 200 ao . The binding energies of the levels are illustrated by 
dashed horizontal lines. 
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